function [ u ] = createPulse( dimR,dimZ, a0, dr, w0,dt,pulseDuration )
    u(dimR,dimZ)=0;
    for r = 1:dimR
        for t=1:dimZ  
        radius = (r-dimR/2)*dr;
        tau = (t-dimZ/2)*dt;
        u(r,t) = a0*exp(-radius*radius/w0/w0)*exp (-tau.*tau/(pulseDuration*pulseDuration));
        end
    end


end

